Transition from fractional to Majorana fermions in Rashba nanowires 
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We study hybrid superconducting-semiconducting nanowires in the presence of Rashba spin-orbit 
interaction (SOI) as well as helical magnetic fields. We show that the interplay between them leads 
to a competition of phases with two topological gaps closing and reopening, resulting in unexpected 
reentrance behavior. Besides the topological phase with localized Majorana fermions (MFs) we find 
new phases characterized by fractionally charged fermion (FF) bound states of Jackiw-Rebbi type. 
The system can be fully gapped by the magnetic fields alone, giving rise to FFs that transmute into 
MFs upon turning on superconductivity. We find explicit analytical solutions for MF and FF bound 
states and determine the phase diagram numerically by determining the corresponding Wronskian 
null space. We show by renormalization group arguments that electron-electron interactions enhance 
the Zeeman gaps opened by the fields. 

PACS numbers: 73.63.Nm; 74.45.+C 



Introduction. Majorana fermions^ (MF) in condensed 
matter systems^, interesting from a fundamental point of 
view as well as for potential applications in topological 
quantum computing, have attracted wide interest, both 
in theory — — and experiment One of the most 
promising candidate systems for MFs are semiconduct- 
ing nanowires with Rashba spin-orbit interaction (SOI) 
brought into proximity with a superconductor^^. In such 
hybrid systems a topological phase with a MF at each end 
of the nanowire is predicted to emerge once an applied 
uniform magnetic field exceeds a critical value — — . As 
pointed out recently^, the Rashba SOI in such wires is 
equivalent to a helical Zeeman term, and thus the same 
topological phase with MFs is predicted to occur in hy- 
brid systems in the presence of a helical field but without 
SOI^^. 

Here, we go a decisive step further and address the 
question, what happens when both fields are present, 
an internal Rashba SOI field as well as a helical-or 
more generally-a spatially varying magnetic field. Quite 
remarkably, we discover that due to the interference 
between the two mechanisms the phase diagram be- 
comes surprisingly rich, with reentrance behavior of MFs 
and new phases characterized by fractionally charged 
fermions (FF), analogously to Jackiw-Rebbi fermion 
bound states^. Since the system is fully gapped by the 
magnetic fields at certain Rashba SOI strengths (in the 
absence of superconductivity), these FFs act as precur- 
sors of MFs into which they transmute by turning on 
superconductivity. 

The main part of this work aims at characterizing the 
mentioned phase diagram. For this we find explicit solu- 
tions for the various bound states, which allows us to de- 
rive analytical conditions for the boundaries of the topo- 
logical phases. We also perform an independent numer- 
ical search of the phases and present numerical results 
illustrating them. We show that the phases can be con- 
trolled with experimentally accessible parameters, such 
as the uniform field or the chemical potential. We for- 
mulate the topological criterion as a condition local in 




FIG. 1. Schematics of the hybrid semiconducting- 
superconducting system, consisting of a finite Rashba 
nanowire (green cylinder) on top of an s-wave bulk supercon- 
ductor (red slab) in the presence of a uniform magnetic field 
B (grey arrow) applied along the nanowire in x-direction. Pe- 
riodically arranged nanomagnets (blue bars) produce a spa- 
tially varying magnetic field Bn(x) (grey arrows). We note 
that Bn(x) can also be generated intrinsically e.g. by a helical 
hyperfine field of nuclear spins inside the nanowire^. 



momentum space via the kernel dimension of the Wron- 
skian, which does not require the knowledge of the spec- 
trum in the entire Brillouin zone. We also address in- 
teraction effects and show that they increase all Zeeman 
gaps and thereby the stability of the topological phase. 

Model. We consider a system consisting of a semicon- 
ducting nanowire with Rashba SOI in proximity with an 
s-wave bulk superconductor and in the presence of mag- 
netic fields which contain uniform and spatially varying 
components, see Fig. [TJ The Rashba spin-orbit interac- 
tion is characterized by a SOI vector a pointing along, 
say, the z-axis. The effective continuum Hamiltonian 
for the nanowire is in Nambu representation given by 
Hq = ^ J dx i/j^ {x)1-LoiIj{x) with 



(1) 



where m is the electron mass. Here, ip = 
(^^,^;,^|,^|), and '^a\x), with a =t / |, is the an- 
nihilation (creation) operator for a spin up/down elec- 
tron at position x. The Pauli matrix {r]i) acts in 
the spin (electron-hole) space. The spectrum of Ho con- 
sists of four parabolas centered at the Rashba momentum 
±/cso = -^ma/h^^ see Fig. [21 The chemical potential /i is 
chosen to be zero at the crossing of the Rashba branches 
at A; = 0. 
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FIG. 2. The spectrum of Rashba nanowire consists of parabo- 
las shifted by ±kso- sohd and dashed hnes correspond to the 
electron and hole spectrum, respectively. The outer circles 
(dotted) mark the exterior branches and the inner circle (dot- 
ted) marks the interior branches. A spatially varying mag- 
netic field with period 4kso (dark green arrow) couples the 
exterior branches. Similarly, a uniform magnetic field couples 
the interior branches ai k — (not indicated). A spatially 
varying magnetic field with period 2kso (light green arrow) 
mixes exterior with interior branches. 

The uniform (B) and spatially varying (B^ ) magnetic 
fields lead to the Zeeman term, 

n, = g^B[B + Bn{x)] ■ a 7^3/2, (2) 

where g is the g-factor and (Ib the Bohr magneton. 
The proximity-induced superconductivity couples states 
of opposite momenta and spins and is described by 1-Ls = 
As (72 7^2, where the effective pairing amplitude Ag can be 
assumed to be non-negative. 

From now on, we assume that the SOI energy mo? jY? 
is the largest energy scale at the Fermi level in the prob- 
lem. In this strong SOI regime, we can treat the 5-fields 
and A5 as small perturbations. This allows us to lin- 
earize the full Hamiltonian Ho + + around = 
(referred to as interior branches) and k = ±2/cso (referred 
to as exterior branches), see Fig. [2l This entails that we 
can use the ansatz 

^(x) = + + L^e-^^^^°^ + i?;e^^^^°^, (3) 

where the right mover Ra{x) and the left mover L^{x) are 
slowly- varying fields. For a uniform magnetic field alone 
(chosen along the x-axis) the full Hamiltonian becomes 
H =\ j dx ^i){x)^l-Lij{x) with 

U = -ihvFcrsrsdx + A^Vs^^i (1 + ^3)/2 + Asr]2(J2, (4) 

where the Pauli matrix acts in the interior-exterior 
branch space, and = (i?^, L^, L|, L^, L|, i?|). 
The Fermi velocity is given by vp = a/h and the Zee- 
man energy by = gfiBB/2. Next, we include the 
spatially varying magnetic field, and assume that it has 
a substantial Fourier component either at Akso (case I) or 
at 2kso (case II), leading to additional couplings between 
all branches of the spectrum, see Fig. [2j We treat now 
the two cases in turn and will see that the interplay of 
Rashba and magnetic fields leads to a surprisingly reach 
diagram of topological phases. 

Case I. Here, we consider Bn{x) with period Akgo 
and perpendicular to a. For a field with oscillating 



amplitude only, we consider two geometries, Bn,x = 
xBn cos{4:ksoX + 0) and Bn^y = yBn sm{4:ksoX + ^), with 
arbitrary phase shift 0, while for a helical field we con- 
sider a field with anticlockwise rotation, (Bn,a; + Bn,2y)/2. 
(clockwise rotation does not lead to coupling). We note 
that Bn{x) can also be generated intrinsically e.g. by 
the hyperfine field of ordered nuclear spins inside the 
nanowire^^. All geometries lead to identical results: they 
affect only the exterior branches (see Fig. [2j) and the 
corresponding Hamiltonian remains block-diagonal in r- 
space. The full Hamiltonian becomes H^kso = ^ + 
An(cri cos6>-cr2 siu 6^)7^3 (1 -r3)/2, where A^ = gjiBBn/^- 
The spectrum for the exterior (/ = e) and interior (/ = i) 
branches is given by 

= {hvphf + A^ + A? + 

I (5) 

where Ag/^ = A^/^. We note the equivalence of effects 
of a uniform field on the interior branches and of a peri- 
odic field on the exterior branches. The spectrum is fully 
gapped except for two special cases, A^^^ = + /i^. 
This suggests that there will be transitions between dif- 
ferent non- trivial phases. 

We identify these phases by the presence or absence of 
bound states inside the gap. For this it is most conve- 
nient to study the Wronskian corresponding to the four 
decaying fundamental solutions^^. Here, we consider a 
semi-infinite nanowire, with boundary at x = 0, and 
assume that all decay lengths will be shorter than the 
system length. For fixed parameters (including the en- 
ergy E)^ we find the four decaying eigenstates of H4kso 
for the left and right movers. Using Eq. (|3|), we express 
them in the basis of the original fermionic fields ?/^, lead- 
ing to four four-spinor solutions with j = 1, 4, and 
construct a 4 x 4 Wronskian matrix Wij{x) = [<l>-^(x)]i. 
The dimension d of the null space of W{0) determines 
the system phase: d = corresponds to a phase with 
no bound states (trivial phase), (i=lat£^ = Otoa 
phase with one single Majorana fermion (MF) (topolog- 
ical phase), (i = 2at£^ = 0or(i=lat£^7^0to 
a phase with one localized fermion of fractional charge 
(FF, see below) (fermion phase). We refer to the triv- 
ial and fermion phases as non-topological. Finally, the 
knowledge of the W{0) null space allows us to construct 
the bound state wave functions, expressed in terms of lin- 
early dependent combinations of fulfilling the Dirich- 
let boundary condition at x = 0. In App. [A] we list the 
analytical solutions for the MF bound states, from which 
we see expHcitly that these solutions are robust against 
any parameter variations (topologically stable) as long as 
the topological gap remains open. 

For case I, we find that the system is in the topological 
phase if one of the following inequalities is satisfied, 

(IA) JR^Al^ < A, < 3?v/Al^ (6) 

(IB) 3?v/Af^ < A, < KVA^ - (7) 



with the corresponding MF wave functions given in 
App. |X1 Case lA goes into IB upon interchange ^ 
A^. As anticipated after Eq. (|5]), the boundaries of the 
topological phase correspond to the system being gapless. 
In the absence of B^, there is only one topological gap, 
which arises from the interior branches^'^. In this case, 
only condition lA can be satisfied, and a MF emerges 
when the uniform 5-field exceeds a critical value. How- 
ever, in the presence of B^, the exterior gap is also topo- 
logical. As shown in Fig.[3l the interplay between the two 
gaps leads to a rich phase diagram with reentrance be- 
haviour. For instance, if \Ai\ > and A^ = 0, the sys- 
tem is in the non-topological phase but still fully gapped 
by the magnetic fields. With increasing A^, first the 
exterior (interior) gap closes and reopens, bringing the 
system into the topological phase. Then, upon further 
increase of Ag, the interior (exterior) gap closes and re- 
opens, bringing the system back into the non-topological 
phase. 

We note that case IB allows the presence of a MF in 
weaker uniform magnetic fields, see Fig. [3l If the nano- 
magnets generating B^ can be arranged such that the 
field penetration into the bulk-superconductor is mini- 
mized, as illustrated in Fig. (U much stronger oscillating 
than uniform fields can be applied, opening up the pos- 
sibility to generate MFs in systems with small ^-factors. 

The system is in the fermion phase, if 



and 



(8) 



where the phases 4>i 
^^i/x) /A 



A?, 



are defined by e*'^^''^ = 
z,n M^+^Mj/^2,n- The corresponding wave 
functions are listed in App. [Al In this regime, two MFs 
(both localized at x = 0) fuse to one fermion bound state. 
Such bound state fermions are known to have fractional 
charge as first discovered in the Jackiw-Rebbi 

model^^^^^. If neither of the inequalities ([6])-(|8]) is satis- 
fied, the system is in the trivial phase without any bound 
state at zero energy. 

Detuning from the conditions in Eq. (j8]) , the two zero- 
energy solutions are usually split, becoming a fermion- 
antifermion pair at energies Importantly, FFs do not 
require the presence of superconductivity. For example, 
if A5 = and /i = 0, the two bound states have energy 



Eff = ± 



A. An sin I' 



VA2 + A2 -2A,An cos^' 



(9) 



We note that the splitting vanishes at 6> = nvr, n integer, 
due to the chiral symmetry of H^kso these special val- 
ues^S. In contrast, the MF remains at zero energy for all 
values of ^, which is a direct manifestation of the stability 
of the bound state within the topological phase (despite 
the fact that the MF wave function depends on 6>, see 
App. [A]). 

To determine the full phase diagram we have per- 
formed a systematic numerical search for all bound state 
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FIG. 3. Phase diagram for case I supporting three phases: the 
trivial phase with no bound states (white area), the topolog- 
ical phase with one MF (green area), and the fermion phase 
with two FFs (colored area). The color in the fermion phase 
encodes the ratio of the fermion energy to the system gap, 
which varies from zero (bright yellow) to one (black). Note 
that the central region of the topological phase corresponds 
to case lA, while the four corner regions to case IB. Here, (a) 
An/ As = 3, 6> = 0, while (b) An/ As = 3, = 7r/4. Note 
that Az < and corresponds to Az > and -\- tt. By 
comparing (a) with (b) and by calculating the dependence on 
[(c) /i/As = 0.4, Az/ As =2], we note that the topological 
phase is insensitive to in contast to the fermion phase. The 
position of a zero-energy FF depends not only on and fi [see 
(a)-(c)] but also on Az and An [(d) /i/As = 1.5, 6 = 7r/4] in 
agreement with Eq. 



solutions with energies inside the gap and determined the 
null space of the Wronskian. The results are plotted in 
Fig. [3l The bright yellow lines inside the colored area in 
Fig. [3] correspond to zero-energy FFs satisfying Eq. (j8j). 
At the point where the lines touch the topological phase 
(shown in green), the gap closes and reopens and the zero- 
energy FF transmutes into a MF. In the fermion phase 
away from the zero-energy line the two solutions split 
(the bigger the splitting the darker the color), until they 
finally reach the gap (black boundaries) and disappear. 
The fermion phase exists only for certain values of the 
phase shift 6>, in contrast to the topological phase, which, 
again, is not sensitive to ^, see Fig. [3t. Moreover, the 
fermion phase is also sensitive to the relative orientation 
of B and B^, see Fig. [3li. In the same panel, we see that 
outwards regions are more suitable for a fractional charge 
observation than the central region, where energies of the 
bound states are very close to the gap edge. 

Case Ila. We now shortly comment on two addititonal 
geometries with a spatially varying magnetic field with a 
period 2kso' To keep the following discussion concise, we 
set /i = and state the results for £^ = bound states 
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only. We begin with a field perpendicular to the SOI 
vector OL and given by B^^^^ = xB^ cos{2ksox) or B^^^ = 
yBn sm{2ksox) (an oscillating field) or (Bn,a; + Bn,2y)/2 (a 
helical field). Such a field mixes the exterior and interior 
branches, see Fig. O The corresponding Hamiltonian is 
^2kso ~ ^ + AnT^sCTiTi. The spectrum is given by 

El = {hvFkf + [Al + Al/4] + (A, ± A,/2)2 (10) 
± 2y^{hvFk)^Al + [Al + A2/4](A, ± A,/2)2. 

Repeating the procedure used above, we derive the 
condition for the topological phase as 



A,>A,|1-A2/a2 



(11) 



Again, the phase boundary to the topological phase cor- 
responds to the parameters at which the system is gap- 
less, i.e. Az = As |l - A^/A^|. We note that in the 
presence of a spatially periodic magnetic field MFs may 
emerge at substantially weaker uniform magnetic fields. 
The fermion phase occurs if 



I An I > A, and A, < A, (A^/A^ - l) 



(12) 



The rest of the parameter space corresponds to the non- 
topological phase. The corresponding wave functions for 
MFs and FFs are given in App. [Al 

Case lib. Finally, we comment on an oscillating field 
aligned with the SOI vector a and given by B^ = 
zBnCos{2ksoX + 6). This field couples the interior and 
exterior branches (see Fig. [2j). The corresponding Hamil- 



tonian is 1-L 
spectrum given by 



2,ks 



H + risasTi cos 9 — T2 sin 9, with the 



Ej = {hvpkf 



(13) 



The topological phase is determined, again, by Eq. (pT]) 
and the corresponding MF wave functions are given in 
App. [Al We note once more that MFs can be observed 
in weaker uniform magnetic fields. Interestingly, in this 
configuration the fermion phase is absent, demonstrating 
the sensitivity of the FFs to the 5-field orientation. 



Interactions. Electron interactions play an important 
role in one-dimensional systems^^ and, in particular, for 
MFs^i^. E.g., the interior gap opened by a uniform 
magnetic field is strongly enhanced by interactions — ; 
and so we expect the same renormalization to occur here 
for both gaps. This is indeed the case, as we show next. 
For this, we perform a renormalization group analysis for 
both the uniform and the periodic field. Following Ref^^ 
we arrive at the effective Hamiltonian H = J ^Hi{x) 
in terms of conjugate boson fields, 0e,i and O^^i^ with 



Hi = v[{\/^iy + (V^O ] + ^ cos(2VK(/)0, (14) 



where we have suppressed quadratic off-diagonal terms 
being less relevant compared to the cosine terms. The 
index / = e/i denotes the exterior/interior branch, a 
the lattice constant, K'^ = {vs/Kg + VcKc)/{vc/Kc + 
VsKs) and v = ^J {vc/Kc + VsKs) {vs/Ks + VcKc)/2. 
Here, Vc^s are the charge (c) and spin (s) velocities 
and Kc^s the corresponding Luttinger liquid parame- 
ters^^. The gaps A^ opened by magnetic fields are 
renormalized upwards by interactions and given by 
Al = Ai{7ThvF/Aia)^^-^y^^-^\ For GaAs (InAs) 



nanowirei 



,21.24 



we estimate an increase by about a fac- 
tor of 2 (4). The enhanced Zeeman gaps allows the use, 
again, of materials with lower ^-factors. 

Conclusions. The interplay between spatially varying 
magnetic fields and Rashba SOI in a hybrid nanowire 
system leads to a rich phase diagram with reentrance 
behavior and with fractionally charged fermions that get 
transmuted into Major ana fermions at the reopening of 
the topological gap. 
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Appendix A: Analytical solutions 

Here we present the analytical solutions of the Ma- 
jorana fermion (MF) and fractional fermion (FF) wave 
functions for different regimes discussed in the main 
text. Every MF wave function has the four-spinor form 
^mf{x) = (/(x), ^(x), /*(x), ^*(x)) with normalization 
condition J dx |^mf(^)P = 2 (below we omit the nor- 
malization factors). For zero energy £^ = 0, the func- 
tions f{x) and g{x) are listed in Table HI together with 
the regime of validity corresponding to the cases I and II 
defined in the main text. 

We use the notations e*"^^ = (y^A^ - /i^ + ifi)/ /S.z, 
e'^- = (v/A^^ + ^/i)/An, and (/)± = (^ ± - (/>.)/2. 
In row Ila-FF of Table [H only one MF wave function 
is given, a second one needs to be added from row Ila- 
MF, depending on whether (A^ - A^/2)2 > A^ or (A^ - 
A./2)^ < A^. 

If A5 = and /i = 0, the fractional fermions have ener- 
gies given be Eq. [9] and the corresponding wave function 
is given by 

= (^^t^ = (^e^Q ^^'^iK.^-^/e. _ e--/«Q , (Al) 

where cosc/)^ = Eff/^z^ = a/\/A^ — Epp, and 
= a/y^A^ — Epp. We note that this fermion of 
non-negative energy has an anti-fermion partner of non- 
positive energy. 
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TABLE I. Wave functions of MPs and FFs for different regimes at zero energy. 



(lA) 
MF 



-^2<A, <7Ai-/i2, \A,,n\>M 

^['^ = a/(A. - VA^-M^), e^^^ = c./(A. + VA^-M^), ^^^^ = a/(yAF^ - A.) 



< A. < VAi-/i2, |An| = ImI, |A,| > 

/ = e-^^-/-+^2fc,ox^_ .g^0./2 _^ (2x/i/a)e-^('/'+"/^) sin[(6> - 0.)/2 + 7r/4]) + ie'^/^'^ +^^^/^ 
g ^ ^-A,x/a-^2k,ox^_^-^c^./2 ^ (2x/i/a) e^^^/^+^'Z^^ sin[(6' - 0^)/2 + 7r/4]) + e— 



< A, < VAi-/i2, |An| < |A,| > = ^/i^ - Al/a 

f = ig*= ^e--/4^^+^^-/2 _ ^g-A,./a+^(2fc,o.+0./2) (^cos(A:^x) + sin(A:^x) [z/i + Ane'^'"^^^] / ^/i^ - A^^ 



(IB) 
MF 



VAi-/i2<A, < |A,,^|>|m| 

/ = -i^* = -ie-^/^i^^ -^^^/^ cos 0- + sin[(^ - 0. + 0.)/2] + ie-^/4'^+^(2^^«^+^-+^^/2) ( 

) = a/(A, - VAi-/i2), = a/(A, + ^AF^), ^i'^ = a/{y^A^^ - A.) 



A. < VAT^, a. = |An| > 

^ ^ g-A,x/a |^•g^(0-0.)/2 _ (2x/,/a)e^"/^ sin [(0n - 6')/2 - 7r/4]) 

^ = e-^^^/- |^_g-.(^-0.)/2 ^ (2^^/c,)e-/4 sin [(0. - ^)/2 - 7t/4]) + e-/4^^+^(4fc...-.+0.)/2 



A, < VA^-M^, |A,| < |A^| > ImI; = ^fi^-Al/a 

f = -ig* = g^(^-0r.)/2 j^.g-x/4^)+z2fc,.x _ ^g-A.x/a (^cos(A:.nx) + sin(A:.nx) [a^c^^^--'^ - i/i] /y^/iS - A^)] 



(I) 

FF 



6* = TT + 0n + 0., < A, < min{3f^7A2 - ^f^^Al"^^}, ^i^i. = 



(Ila) 
MF 



A, > A, |1 - A^/AJI and (A, - A,/2)2 > A^; = (A./2 ± ^(A. - A,/2)2 - A2 

f = ig*= ie^-f^^ox (^e— + e-'/^+A^ - ie-^/^+An - ie""/^- A, A = a/^- - A, 



A, > A, |l - A^/A^l and (A, - A,/2)^ < Al; cos0 = (A,/2 - A3)/An, km = y/^n - (A. - A,/2)Va 
f = ig*= [ie^^^-^ cos(/c,nX + 0/2) - i cos(A:,nX - 0/2)] g-^^'/^" 



(Ila) 
FF 



|An| > A., A. < A. (A^/A^ - 1), and (A, + A,/2)2 > A^ ; = a/ (a./2 ± ^(A. + A./2)2 - A^ 



|An| > A., A. < A. (A^/Aj - 1) , and (A, + A,/2f < Al; cos = (A,/2 + A,)/An, km = ^A^ - (A. + A./2)2 
/ = -ig* = [cos{kmX - 0/2) - e^^^^^-" cos{kmX + 0/2)] e-^-"/^— 



(lib) 
MF 



A, > A, |A^/Aj - l| 

+i[A,/2 + VA2 + Ai/4](e^(''^^°"-')(e-"/^2 ^08 6* + ie-"/^^ sin (9) - e'^/^^) 
= -iAn[e-^/^2 cos6> + sin 6* - e-^/^3-^(2fc.oa.-0)] 



+ [A,/2 + ^Al + Ai/4](e-^(2'^^°"-^He-"/^2 cos6> - ie-"/^^ sin 6*) - e'^/^^) 
6 = a/ (a./2 + a. - VA2 +Ai/4) , 6 = «/ (a./2 - A. + y^A^ + A2/4) , ^4 = «/ (a. + ^Al + Ai/4 - A./2) 



